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ORIGINAL ARTICLE 

Kinetic Interpretation of the Importance of OATP1B3 and 
MRP2 in Docetaxel-lnduced Hematopoietic Toxicity 

A Yamada 1 , K Maeda 1 , K KiyotanP, T Mushiroda 2 , Y Nakamura 43 and Y Sugiyama 5 

Neutropenia is a lethal dose-limiting toxicity of docetaxel. Our previous report indicated that the prevalence of severe docetaxel- 
induced neutropenia is significantly associated with genetic polymorphisms in solute carrier organic anion transporter 1B3 
(SLC01B3) (encoding organic anion-transporting polypeptide 1 B3 (OATP1 B3)) and ATP-binding cassette subfamily C2 (ABCC2) 
(encoding multidrug-resistant-associated protein 2 (MRP2)). Therefore, we investigated their significance in docetaxel-induced 
neutropenia. In vitro experiments suggested their possible involvement in the hepatic uptake of docetaxel and its efflux 
from bone marrow cells. To further characterize a quantitative impact of OATP1B3 and MRP2 on neutropenia, we used an in 
siiico simulation of the neutrophil count in docetaxel-treated subjects with functional changes in OATP1B3 and MRP2 in a 
pharmacokinetic/pharmacodynamic model. The clinically reported odds ratios for docetaxel-induced neutropenia risk were 
explained by the decreased function of OATP1B3 and MRP2 to 41 and 32%, respectively. These results suggest that reduced 
activities of OATP1B3 and MRP2 associated with systemic exposure and local accumulation in bone marrow cells, respectively, 
account for the docetaxel-induced neutropenia observed clinically. 
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INTRODUCTION 

Docetaxel exerts an antitumor activity by stabilizing and 
inhibiting the depolymerization of tubulin, 1 and is used as a 
first-line therapy in the treatment of several kinds of cancers. 
Neutropenia is a dose-limiting toxicity of docetaxel, and it 
sometimes restricts the clinical use. However, the determin- 
ing factors for docetaxel-induced hematopoietic toxicity are 
not understood well. 

When infused intravenously in humans, docetaxel is 
metabolized mainly by cytochrome P450 (CYP) 3A in the 
liver, and only a small fraction of the parent docetaxel is 
excreted into the bile and urine. 2 Because the metabolites of 
docetaxel have less biological activity than the parent form, 3 
systemic exposure to docetaxel is thought to be an important 
causal factor of docetaxel-induced toxicity. Previous reports 
suggested that the area under the plasma concentration- 
time curve of docetaxel correlated with a decrease in neu- 
trophil count. 45 

Kiyotani etal. 6 reported that the risk of severe neutropenia 
induced by docetaxel chemotherapy was significantly asso- 
ciated with genetic polymorphisms of solute carrier organic 
anion transporter 1B3 (SLCOW3), encoding organic anion- 
transporting polypeptide 1 B3 (OATP1 B3), and of ATP-binding 
cassette subfamily C2 (ABCC2), encoding multidrug-resis- 
tant-associated protein 2 (MRP2). It was reported earlier 
that docetaxel was a substrate of OATP1B3 and MRP2. 7-10 
OATP1 B3 is expressed almost exclusively on the basolateral 
side of human hepatocytes and mediates the hepatic uptake 
of a wide variety of drugs such as telmisartan. 11-13 Although 
docetaxel is reported to be a substrate of OATP1B3, the 



relative contribution of OATP1B3 to hepatic uptake of 
docetaxel remains to be clarified. MRP2 is expressed on the 
bile canalicular membrane and is responsible for the biliary 
excretion of various organic anions including glucuronide and 
glutathione conjugates. 14-16 However, because docetaxel is 
eliminated predominantly by cytochrome P450 3A-mediated 
metabolism, it is unlikely that hepatic MRP2 controls its sys- 
temic exposure. Instead, MRP2 in hematopoietic cells or their 
precursor cells is likely to suppress the entry of docetaxel into 
cells and protects them against its cytotoxicity. 

The pharmacokinetic/pharmacodynamic (PK/PD) model is 
a powerful tool for predicting both drug efficacy and toxicity. 
A population PK/PD model was built to predict the time pro- 
file of blood neutrophil count after treatment with anticancer 
drugs including docetaxel, and this model could explain well 
the quantitative risk of neutropenia. 17 Because the mutations 
in SLCOW3 and ABCC2 reported previously are not located 
in their coding regions, it is difficult to estimate the functional 
alterations in these transporters from in vitro experimental data. 
Alternatively, we decided to use the reported PK/PD model 
describing the drug-induced neutropenia for that purpose. 

In this study, we tried to investigate the possible roles of 
OATP1B3 and MRP2 in the systemic pharmacokinetics and 
local exposure of docetaxel in bone marrow cells, which finally 
links to the hematopoietic toxicity. Contribution of OATP1B3 
to the overall uptake of docetaxel in human hepatocytes and 
the protective role of MRP2 in the cytotoxicity induced by 
docetaxel were investigated by in vitro experiments. Then, 
the impact of functional changes in OATP1 B3 and MRP2 on 
docetaxel-induced toxicity was also estimated quantitatively 
using Monte Carlo simulation approach. 
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RESULTS 

Uptake of docetaxel in hepatic uptake transporter 
expression systems 

Among hepatic uptake transporters in humans (OATP1B1, 
OATP1B3, OATP2B1 , OAT2, Na + -taurocholate cotransport- 
ing polypeptide (NTCP), and organic cation transporter 
(OCT)1), docetaxel is taken up significantly only into 
OATP1B3-expressing HEK293 cells compared with con- 
trol cells (Figure 1a). Transport activity of each expression 
system was confirmed as the uptake of typical substrates, 
estrone-3-sulfate (E^) for OATP1B1 and OATP2B1 (49 
and 14 ul/0.5min/mg protein), cholecystokinin octapep- 
tide (CCK-8) for OATP1B3 (35 ul/5min/mg protein), xan- 
thine for OAT2 (18 ul/15min/mg protein), taurocholate for 
NTCP (14 ul/2min/mg protein), and tetraethylammonium 
for OCT1 (90 ul/10min/mg protein). The saturation kinet- 
ics of OATP1B3-mediated uptake of docetaxel was evalu- 
ated by its transport for 5min because the time-dependent 
linear uptake was maintained (data not shown). The pro- 
tein-unbound fraction of docetaxel (0.001-100 umol/l) in 
the transport buffer with 3% human serum albumin was 
measured to estimate the transport clearance with regard 
to the unbound concentration of docetaxel. The protein- 
unbound fraction was constant within the range of 0.001-1 
umol/l, but increased gradually in the concentration range 
of 3-100 umol/l, suggesting a saturation of protein bind- 
ing (Supplementary Figure S1). The uptake clearance 



was calculated with regard to the unbound concentration 
and is shown as Eadie-Hofstee plots (Figure 1b). The 
concentration-dependent uptake could be explained by 
one saturable component with a K m value of 0.325 ±0.064 
umol/l (mean ± computer-calculated SD). 

Inhibitory effects of E^S and estradiol-1 7p-glucuronide 
on docetaxel uptake 

The in vitro uptake study with OATP1B1- and OATP1B3- 
expressing HEK293 cells confirmed that 2 mmol/l E A S 
with 3% human serum albumin specifically inhibited 
OATP1B1 -mediated transport but not OATP1 B3-mediated 
transport (Figure 1c,d). By contrast, 2 mmol/l estradiol- 
17p-glucuronide (E 2 17pG) inhibited both OATP1B1- and 
OATP1B3-mediated transport. The uptake of docetaxel in 
human hepatocytes was inhibited by 100 umol/l unlabeled 
docetaxel and 2 mmol/l E 2 17pG, but was not inhibited by 2 
mmol/l E^S (Figure 1e). 

Cytotoxicity and intracellular accumulation of docetaxel 
in MRP2-expressing MDCKII cells 

The LC 50 values in MRP2-expressing and vector-trans- 
fected MDCKII cells were calculated to be 12.0±5.4 and 
3.27 ±0.91 nmol/l, respectively (mean ± computer-calcu- 
lated SD; Figure 2a). The accumulation of docetaxel for 
24 h in these cell lines was also measured. The distribution 
volume in MRP2-expressing cells (271 ±9 ul/mg protein) 




Figure 1 In vitro uptake of docetaxel into HEK293 cells expressing solute carrier transporters expressed in human hepatocytes and into human 
hepatocytes. (a) The uptake of docetaxel (1 umol/l) into OATP-, OAT-, OCT1-, and NTCP-expressing HEK293 cells and vector-transfected 
control cells (vector) was measured at 5min. (b) Saturation of the docetaxel uptake in OATP1 B3-expressing HEK293 cells was also investigated. 
The solid line is a fitted curve calculated by nonlinear regression analysis based on Eq. 1 , as described in the Methods, (c) Concentration (0.1 
umol/l and 2 mmol/l)-dependent uptake of estrone-3-sulfate (E^) in OATP1 B1 -expressing HEK293 cells, (d) inhibitory effects of 2 mmol/l 
estradiol-1 7p-glucuronide (E 2 17pG) and E^ on the uptake of cholecystokinin octapeptide (CCK-8) in OATP1 B3-expressing HEK293 cells, 
and (e) docetaxel uptake into cryopreserved human hepatocytes in the presence of E^S and E 2 17pG were investigated. These uptake assays 
were carried out in the presence of 3% human serum albumin. Each bar represents the mean ± SE (n = 3). *P< 0.05, **P< 0.01 , ***p< 0.001 . 
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Figure 2 Effect of exogenous expression of MRP2 on the cytotoxicity 
(a) and intracellular accumulation (b) of docetaxel. (a) Vector- 
transfected MDCKII cells (closed circles) or MRP2-expressing 
MDCKII cells (open circles) were incubated for 96 h with various 
concentrations (0.01-100 nmol/l) of docetaxel. Cell survival was 
quantified by the MTT assay as described in the Methods, (b) Cell 
lines were incubated for 24 h with 1 nmol/l of docetaxel, and the 
volume uptake of docetaxel into these cells was determined. Each 
point and bar represents the mean ± SE (n = 3). *P< 0.05, **p< 
0.01 ,***P< 0.001. 

was significantly smaller than that in vector-transfected cells 
(665 ±20 ul/mg protein) (Figure 2b). 

Colony-forming assay 

A colony-forming assay was performed to investigate the inhibi- 
tory effect of docetaxel on the proliferation and maturation of 
bone marrow cells induced by granulocyte colony-stimulating 
factor. Docetaxel decreased the number of colonies in a con- 
centration-dependent manner. The inhibitory effect of docetaxel 
on colony formation was significantly more potent in bone mar- 
row cells derived from the MRP2-deficient rats (Eisai hyper- 
bilirubinemic rats), compared with cells from normal Sprague 
Dawley rats (Figure 3a). The inhibition potency of docetaxel 
was also significantly increased in the presence of 50 umol/l of 
the MRP inhibitor, MK571 (Figure 3b). 

Estimation of EC 50 values for subjects with no targeted 
mutations in SLC01B3 and ABCC2 for an in silico simu- 
lation using the PK/PD model 

We initially performed an in silico simulation by using 
a reported PK/PD model describing the time course of 
docetaxel pharmacokinetics and neutrophil count (Figure 5; 
Supplementary Material online). 1718 However, in subjects 
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Figure 3 Inhibitory effect of docetaxel on granulocyte colony- 
stimulating factor (G-CSF)-induced colony formation, (a) Bone 
marrow cells obtained from Sprague Dawley rats (open circles) and 
Eisai hyperbilirubinemic rats (closed circles) were incubated for 6 
days with various concentrations (0-5 ng/ml) of docetaxel. (b) Bone 
marrow cells derived from Sprague Dawley rats were also incubated 
for 6 days with various concentrations (0-5 ng/ml) of docetaxel in 
the presence (closed circles) or absence (open circles) of 50 umol/l 
MK571 . The number of colonies induced by G-CSF was counted. 
Each point represents the mean ± SE (n = 9).*P < 0.05, **P < 0.01 , 
***P< 0.001. 

with wild-type (wt) SLCOW3 and ABCC2 alleles, the ratio of 
the number of patients with grade 3-4 neutropenia to those 
without neutropenia calculated from an in silico simulation 
with reported PK/PD parameters (1.71) largely differed from 
the clinically observed data (0.0968) (data not shown). We 
considered that this difference might reflect the difference 
in the EC 50 value in this model because the very frequent 
ABCC2 genotype, which can affect the transport function of 
MRP2 and subsequently the EC 50 value, was not taken into 
consideration in the previous parameter estimation. 17 There- 
fore, we re-evaluated the EC Kn value that best fitted our clini- 
cal data. Setting the EC 50 values to 10, 20, 22, 24, 26, 28, 
or 30 umol/l gave the ratios of the number of patients with 
grade 3-4 neutropenia to those without neutropenia of 1 .09, 
0.178, 0.146, 0.146, 0.106, 0.0880, and 0.0864, respectively 
(Figure 4a). This estimation gave an optimum EC 50 value that 
would best account for our clinical data of about 28 umol/l. 

Simulated odds ratios in subjects with mutations in 
SLC01B3 and ABCC2 

In the heterozygotes of SLCOW3 rs1 1045585 (A>G), the 
intrinsic hepatic clearance was expected to decrease to 
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Figure 4 Estimation of EC 50 values in the control group and the 
functional changes to MRP2, which best explained the clinical 
observations reported previously. 6 (a) Ratio of the number of 
subjects with severe neutropenia (grade 3 or 4) to those without any 
hematopoietic toxicity (non-ADR) was calculated with EC 50 values 
of 10, 20, 22, 24, 26, 28, or 30 umol/l, respectively. The broken 
line shows the clinically reported value (0.0968). (b) The impact of 
MRP2 function on the OR for the expression of severe (grade 3 or 
4) neutropenia was simulated. The upper, middle, and lower broken 
lines represent the ORs for OATP1B3 (wt/wt)/MRP2 (mut/mut) 
(1 1 .8), OATP1 B3 (wt/mut)/MRP2 (wt/mut) (9.04), and OATP1 B3 (wt/ 
wt)/MRP2 (wt/mut) (3.72), respectively, as observed clinically. ADR, 
adverse drug reaction; OR, odds ratio. 

70.5% because this was estimated from the difference in the 
intrinsic hepatic clearance of telmisartan between homozy- 
gotes of wt SLCOW3 alleles (wt/wt) and heterozygotes of 
the mutated (mut) SLCOW3 allele (wt/mut) in our previous 
pharmacogenetic study. 19 On the basis of this information, 
the hepatic clearances of docetaxel were estimated to be 
36.7, 29.5, and 20.0 l/h in subjects with wt/wt, wt/mut, and 
homozygotes of mutated alleles (mut/mut) of SLCOW3, 
respectively, assuming that the change in the intrinsic clear- 
ance of telmisartan was the same as that of docetaxel. 

Regarding the functional change of MRP2 caused by the 
ABCC2 rs1 2762549 (G>C) mutation, we estimated the odds 
ratio (OR) in subjects with mutated alleles of ABCC2 by 
changing the EC 50 value when the in vivo transport activity 
of MRP2 from mutated ABCC2 was assumed to decrease to 
10-90% (Figure 4b). By minimizing the sum of the square 
of differences normalized by observed ORs between clini- 
cally observed ORs in subjects with SLCOW3 and ABCC2 
diplotypes of wt/wt and wt/mut (OR wt/wt+wt/mut ; 3.72), wt/wt and 
mut/mut (OR wt/wt+mut/mut ; 11.8), and wt/mut and wt/mut (OR wt/ 
mut+wt/mut' 9 - 04 ) ancl calculated 0 Rs with changing EC 50 value, 



Table 1 Comparison of the number of patients who received docetaxel and 
the ORs for docetaxel -induced neutropenia between clinically observed data 
and simulated data (each genotype: n = 500) when the mean EC 50 value was 
set at 28 umol/l. 
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Non-ADR, nonadverse drug reaction; OR, odds ratio. 

we conclude that the decrease in the MRP2 function to 32% 
would best explain our clinical data (Table 1). 

DISCUSSION 

Although docetaxel is used in the treatment of several kinds 
of cancers, neutropenia sometimes impedes its continuous 
clinical use. The severity of neutropenia is thought to be 
determined by the exposure to docetaxel in the bone mar- 
row, which is dominated by both its systemic concentration 
and local accumulation in hematopoietic cells or their pre- 
cursors. In the present study, we tested our hypothesis that 
hepatic uptake of docetaxel is controlled mainly by OATP1 B3 
and that MRP2 in hematopoietic cells acts as a barrier to the 
entry of docetaxel into hematopoietic cells. 

Among the hepatic uptake transporters, only OATP1B3 
recognizes docetaxel as a substrate, as reported previ- 
ously 78 ; however, its relative contribution to the hepatic 
uptake of docetaxel is unknown. The uptake of docetaxel into 
human hepatocytes was inhibited by E 2 17pG (an inhibitor of 
both OATP1B1 and OATP1B3), but was not inhibited by E^S 
(an OATP1 B1 -selective inhibitor). 11 This suggests that the 
hepatic uptake of docetaxel is mediated mainly by OATP1 B3, 
as is the case for telmisartan. These results support our 
hypothesis that the increased risk of docetaxel-induced neu- 
tropenia by the SLCOW3 rs1 1045585 mutation is caused by 
the decreased expression of OATP1B3 because this muta- 
tion is not located in the exon of SLC01B3. 

Because the fraction of biliary excretion of parent 
docetaxel is thought to be very small, we hypothesized that 
MRP2 might play an important role in bone marrow cells. 
Previous reports showed that MRP2 expression was not 
detected in CD34 _ KSL cells (hematopoietic stem cells), but 
was detected in CD34 + KSL cells (immature hematopoietic 
cells). 20 In human MRP2-expressing cells, the intracellular 
accumulation of docetaxel decreased significantly, and sen- 
sitivity to the cytotoxicity defined based on its medium con- 
centration subsequently decreased. Moreover, bone marrow 
cells from Mrp2-deficient Eisai hyperbilirubinemic rats were 
more sensitive to docetaxel-induced cytotoxicity compared 
with cells from normal Sprague Dawley rats. The presence 
of MK571, an MRP inhibitor, also sensitized bone marrow 
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Figure 5 Schematic diagram of the PK/PD model describing docetaxel-induced neutropenia, (a) A three-compartment model was used 
to describe the pharmacokinetics of docetaxel. (b) The model of Friberg etal. was used to estimate the neutrophil count. 17 (c) Neutrophil 
concentration profile after docetaxel administration in five individuals selected from subjects with SLCOW3 and ABCC2 diplotypes of wt/wt 
and wt/wt. Broken lines show borders between non-ADR, neutropenia grade 1/2, and neutropenia grade 3/4. Non-ADR, nonadverse drug 
reaction; PD, pharmacodynamic; PK, pharmacokinetic. 



cells to docetaxel. These results support our hypothesis that 
MRP2 in bone marrow cells works as a barrier to docetaxel- 
induced hematopoietic cytotoxicity. 

The single-nucleotide polymorphisms in SLCOW3 and 
ABCC2 linked to the increased risk of docetaxel-induced 
neutropenia are located in the noncoding region of each gene 
(intron 1 1 and the 3'-untranslated region, respectively). 6 Thus, 



it is difficult to estimate the quantitative functional changes of 
these transporters from an in vitro study, although we can- 
not neglect the possibility that unknown causal mutations are 
linked to these single-nucleotide polymorphisms. Telmisartan 
is expected to be an in vivo probe substrate for OATP1 B3. 1121 
We demonstrated recently that the intrinsic clearance of 
telmisartan decreased to 70.5% in heterozygotes with the 
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SLCOW3 rs1 1045585 mutation compared with homozy- 
gotes of wt alleles. 19 That result fits our hypothesis, explain- 
ing why the risk of docetaxel-induced hematopoietic toxicity 
is increased by the SLCOW3 mutation. On the other hand, 
because MRP2 is not involved in the biliary excretion of 
docetaxel, decreased expression of MRP2 might occur in 
certain types of bone marrow cells. 

A mathematical PK/PD model for predicting neutropenia 
induced by anticancer drugs including docetaxel has been 
proposed. 17 In this model, a multicompartment PK model 
is connected to the PD model that is used to describe the 
maturation, feedback regulation, and turnover of hemato- 
poietic cells and their precursor cells. The elimination rate 
constant (/C 10 ) of docetaxel is affected directly by a functional 
change in OATP1 B3. By contrast, the EC 50 value was defined 
with regard to the intracellular concentration of docetaxel in 
hematopoietic cells in this model, and so the decreased func- 
tion of MRP2 corresponds to the decrease in EC 50 values. 

To understand the distribution of the clinical outcomes 
in a certain population, a set of parameters was generated 
randomly in the PK/PD model for each subject (a so-called 
"virtual person") while retaining the distribution of parameters 
in one population using Monte Carlo simulation (Figure 5c). 
Kato et al. used this approach to demonstrate that the inter- 
individual variability of cytochrome P450 3A4 function could 
be estimated from the variation in the clinical pharmacoki- 
netics of cytochrome P450 3A4 substrates. 22 In our simu- 
lation, assuming that the variability of each parameter was 
unchanged when the parameter (/C 10 or EC 50 ) reflecting the 
transporter function was changed by genetic polymorphisms 
of SLCOW1 and ABCC2, the model could predict the distribu- 
tion of neutrophil counts in subjects with different genotypes. 

Initially, using the default parameters reported, 17 the cal- 
culated ratio of the number of patients with severe (grade 
3-4) neutropenia to that without any hematopoietic toxicity 
in the control group was 1.71, which differed considerably 
from the reported value of 0.0968 (data not shown). One pos- 
sible reason is the high frequency of the ABCC2 rs1 2762549 
genotype in the Japanese population (0.636). To fit the real 
situation, we first optimized the EC 50 value to adjust it to 
the reported value. When the EC 50 value in subjects with wt 
SLCOW3 and ABCC2 alleles was set at 28 umol/l, the calcu- 
lated ratio was 0.0880, which was close to the reported value. 

The decrease of oral clearance in heterozygotes of SLCOW3 
rs1 1045585 could explain the clinically observed OR for the 
risk of docetaxel-induced neutropenia. A metabolomic study 
revealed that the plasma concentrations of some MRP2 sub- 
strates were increased in carriers of ABCC2 rs12762549, 23 
suggesting a decreased function of MRP2. However, it is dif- 
ficult to determine the extent of its functional change. There- 
fore, we decided to estimate the function in mutated MRP2 by 
comparing the reported ORs for neutropenia with the calcu- 
lated ones in subjects with SLCOW3 and ABCC2 genotypes 
of wt/mut and wt/mut, wt/wt and wt/mut, and wt/wt and mut/mut, 
with numbers of more than 15 in each group. When the MRP2 
function was assumed to decrease to 32% by this mutation, the 
calculated OR was similar to the reported value. Our sensitivity 
analysis showed that even a small change in MRP2 function 
has a large effect on the OR. Previous reports suggested that 
many factors, such as environmental and pathophysiological 



conditions, can change the expression of MRP2 24 and that the 
large interindividual differences in the risk of neutropenia partly 
reflect the differences in the accumulation of docetaxel into cer- 
tain types of bone marrow cells, which are partly determined 
by MRP2. 

In conclusion, we demonstrated that the hepatic uptake of 
docetaxel is mediated mainly by OATP1B3 and that MRP2 
plays an important role as a barrier against the entry of 
docetaxel into hematopoietic cells.The insilicoPK/PD analyses 
showed that the increased risk of docetaxel-induced neutro- 
penia could be explained by functional decreases in OATP1 B3 
and MRP2. Our model could explain the clinically observed 
ORs for docetaxel-induced neutropenia when the function of 
OATP1B3 and MRP2 decreased to 41 and 32%, respectively, 
because of their single-nucleotide polymorphisms. 

METHODS 
Materials 

[ 3 H]-docetaxel (60 Ci/mmol) and [ 3 H]-xanthine (19.4 Ci/ 
mmol) were purchased from American Radiolabeled Chemi- 
cals (St Louis, MO). [ 3 H]-E 2 17pG (45 Ci/mmol), [ 3 H]-E£ (46 
Ci/mmol), [ 3 H]-taurocholic acid (5.0 Ci/mmol), and [ 14 C]- 
tetraethylammonium (3.2 mCi/mmol) were purchased from 
PerkinElmer Life and Analytical Sciences (Boston, MA). [ 3 H]- 
CCK-8 (77 Ci/mmol) was purchased from Amersham Phar- 
macia Biotech (Buckinghamshire, UK). Unlabeled E 2 17pG, 
E^S, and CCK-8 were purchased from Sigma-Aldrich (St 
Louis, MO), and docetaxel was purchased from LC Labora- 
tories (Woburn, MA). All other chemicals were of analytical 
grade and are commercially available. 

Animals 

Male Sprague Dawley rats and Eisai hyperbilirubinemic rats 
(7-8 weeks old) were used in this study. All animals were 
maintained under standard conditions with a reversed dark- 
light cycle and were treated humanely. Food and water were 
available ad libitum. The studies were carried out in accor- 
dance with the guidelines provided by the Institutional Animal 
Care Committee (Graduate School of Pharmaceutical Sci- 
ences, the University of Tokyo, Tokyo, Japan). 

Transport studies using transporter expression systems 

Human OATP-, OAT-, or NTCP-expressing HEK293 cells 
have been constructed previously. 25-27 Human OCT1-trans- 
fected HEK293 cells were constructed as follows. Human 
OCT1 (GenBank Accession No. NM_153187) transcripts 
were amplified and subcloned into pcDNA3.1 vector (Invi- 
trogen, Carlsbad, CA). The expression vector carrying 
human OCT1 cDNA was then transfected into HEK293 
cells using Lipofectamine 2000 (Invitrogen) according to 
the manufacturer's protocol and stably transfected cells 
were selected by adding G418 (0.8mg/ml) (Invitrogen) to 
the culture medium. The transport study was carried out as 
described previously. 27 

Transport study using cryopreserved human hepatocytes 

Cryopreserved human hepatocytes (Lot. 109) were pur- 
chased from BD Gentest (Franklin Lakes, NJ). Immediately 
before the study, the hepatocytes in 1-ml suspension were 
thawed at 37°C and then quickly suspended in thawing 
medium (Celsis In Vitro Technologies, Baltimore, MD). The 
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cells were centrifuged, the medium was replaced with plat- 
ing medium, and the hepatocytes were transferred to col- 
lagen-precoated 24-well plates (Becton Dickinson Labware, 
Bedford, MA) at a cell density of 1 .8x 10 5 viable cells/well in 
0.5 ml of plating medium and then incubated for 24 h. The 
uptake study was performed as described previously. 27 

Kinetic analysis of OATP1 B3-mediated transport 

Docetaxel uptake was expressed as volume (ul/mg protein) 
and given as the amount of radioactivity associated with the 
cells (dpm/mg protein) divided by its initial concentration in the 
incubation medium (dpm/ul). Specific uptake was calculated 
by subtracting uptake in the vector-transfected cells from that 
in the transporter cDNA-transfected cells. The kinetic param- 
eters were calculated using the following equation: 



v = ¥nvrS_ + p 



K m +S 



dif " 



0) 



where v, S, K m , V max , and P dif represent the uptake velocity 
of the substrate (pmol/min/mg protein), substrate concen- 
tration in the medium (umol/l), Michaelis constant (umol/l), 
maximum uptake rate (pmol/min/mg protein), and nonsatu- 
rable uptake clearance (ul/min/mg protein), respectively. Fit- 
ting was performed using a nonlinear least-squares method 
using the MULTI program, 28 and the Damping Gauss-Newton 
algorithm was used for curve fitting. These input data were 
weighted as the reciprocals of the observed values. 

To measure the saturable hepatic uptake clearance in 
human hepatocytes, we first determined the hepatic uptake 
clearance (CL 



(2-0.5 min)' 



(ul/min/mg protein) by calculating the 
slope of the uptake volume (V 6 ) (ul/mg protein) between 0.5 
and 2min (Eq. 2). The saturable component of the hepatic 
uptake clearance (CL hep ) was determined by subtracting 
CL (2 05min) in the presence of 100 umol/l substrate (excess) 
from that in the presence of 10 nmol/l (E^ and CCK-8) or 
10.9 nmol/l (docetaxel) substrate (tracer) (Eq. 3): 



CL 



(2-0.5min) ' 



2-0.5 



CL 



hep 



CL, 



(2-0.5min),tracer 



■CL 



(2-0. 5min), excess 



(2) 



(3) 



MTT viability assays 

Human MRP2-expressing and vector-transfected MDCKII 
cells were seeded in 96-well plates at 1x10 3 cells/well. After 
12h incubation, the medium was replaced with Dulbecco's 
modified Eagle's medium containing docetaxel at concentra- 
tions of 0.01-100 nmol/l. After 96h, the medium was removed 
and the cells were washed with phosphate-buffered saline 
twice. Fresh Dulbecco's modified Eagle's medium (without 
phenol red) containing 100 ul of 3-(4,5-dimethylthiazol-2-yl)- 
2,5-diphenyltetrazolium bromide (MTT) was added to each 
well. After 3h of incubation, the medium was replaced, 100 ul 
of dimethyl sulfoxide was added to each well to dissolve the 
formazan, and the absorbance at 560 nm was measured. The 
LC 50 value was calculated using the following equation: 



Survival(%)(/) = 



Survival(%)(/ = 0) 



1 + 



/ 



(4) 



where / represents the concentration of docetaxel. Fitting 
was performed with a nonlinear least-squares method using 
the MULTI program, 28 and the Damping Gauss-Newton algo- 
rithm was used for curve fitting. 

Accumulation study 

MRP2-expressing and vector-transfected MDCKII cells 
were seeded in 12-well plates at 1.25x10 5 cells/well. 
After 48-h incubation, the medium was replaced with fresh 
medium containing [ 3 H]-docetaxel and the cells were incu- 
bated for 24 h. The cells were then washed three times with 
1 ml of ice-cold phosphate-buffered saline, solubilized with 
500 ul of 0.2 N NaOH and stored overnight at 4 °C. After 
250 ul of 0.4 N HCI was added to solubilized cells, aliquots 
(600 ul) were transferred to scintillation vials. The radioac- 
tivity associated with the cells and the incubation buffer was 
measured with a liquid scintillation counter after the addi- 
tion of 3ml of scintillation fluid (Clear-sol I; Nacalai Tesque, 
Kyoto, Japan) to the scintillation vials. 

Colony-forming assay 

The inhibitory effect of docetaxel on the colony formation of 
bone marrow cells induced by granulocyte colony-stimulat- 
ing factor was performed as described. 29 Bone marrow cells 
were collected from the femurs of rats by flushing dissected 
and cleaned femurs with a-medium (GIBCO BRL, Gaith- 
ersburg, MD) containing 10% fetal bovine serum. The cells 
were collected and then incubated for 2h. The nonadher- 
ent cells were passed through a 40-um cell strainer (Becton 
Dickinson Labware, Franklin Lakes, NJ) to prepare single 
cell suspensions. The methylcellulose colony-forming assay 
was performed in 35-mm dishes. The culture medium con- 
sisted of a-medium, 1.2% 1500-cps methylcellulose, 30% 
fetal bovine serum, 10% bovine serum albumin, 96 umol/l 
2-mercaptoethanol, 0.1 mmol/l 1-thioglycerol, and 0.1% 
Neutrogin (Chugai Pharmaceuticals, Tokyo, Japan). Cells 
were plated in 35-mm dishes at a concentration of 1 x10 6 
cells and incubated for 6 days. Aggregates comprising more 
than 50 cells were counted as colonies and quantified by 
observation with an inverted microscope (Diaphot; Nikon, 
Tokyo, Japan). 

In silico simulation of PK/PD model for predicting neu- 
trophil count after the treatment of docetaxel. 

The reported PK/PD model shown in Figure 5 was used 
to predict the neutrophil count after docetaxel administra- 
tion. 1718 The differential equations used in the PK/PD model 
were as follows: 



5 = /-(^10+^12+^13)xX 1 +/( 21 xX 2 +/(3 1 xX3 (5) 

at 



dt 

dXs 
dt 



'■ ^13 X ^1 _ ^31 X ^3 
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Table 2 Descriptions for all the parameters and mean values and their CV of kinetic parameters in the PK/PD model for the in silico simulation of docetaxel- 
induced neutropenia. 1718 



Parameter 


Description 


Mean value 


© (CV%) 


/(umol/h) 


Infusion rate of docetaxel 


130 




X 


Amount of docetaxel in compartment /'in Figure 5a 


Simulated 


— 


K 10 (h- 1 ) 


Elimination rate constant of docetaxel 


Calculated (Eq. 8) 


— 


K 12 (h- 1 )* 


Distribution rate constant from compartment 1 to 2 


1.07 


— 


K 21 (h- 1 ) a 


Distribution rate constant from compartment 2 to 1 


1.74 


131 


k 13 (h-y 


Distribution rate constant from compartment 1 to 3 


0.0787 


14.7 


K 31 (h- 1 ) a 


Distribution rate constant from compartment 3 to 1 


1.28 


47.7 


CL tot (l/h) a 


Systemic clearance of docetaxel 


36.7 


33.5 


V, (l) a 


Volume of central compartment 


8.31 


56.1 


Prol 


Compartment representing the number of hematopoietic stem cells 
and progenitor cells 


Simulated 





Transit 


Transit compartment describing the maturation of hematopoietic cells 


Simulated 


— 


*p ro , (h- 1 ) 


Proliferation rate constant determining the rate of cell division 


Calculated (Eq. 11) 


— 


Circ (x 10 9 /l) 


Compartment representing the number of neutrophils in the circulat- 
ing blood 


Simulated 


— 


Circ 0 (x 10 9 /l) b 


Baseline neutrophil concentration 


5.05 


42 


(CircyCirc^ 


Feedback regulating the number of circulating neutrophils 


Simulated 




f 


Hill coefficient of negative feedback on cell proliferation 


0.163 


— 


MTT (h) b 


Mean transit time 


89.3 


16 


N 


Number of transit compartments 


3 




^Drug 


Drug-induced effect described with E max model 


Simulated 






Plasma concentration of docetaxel 


Simulated 




MM 


Maturation rate constant representing the time delay 


Calculated (Eq. 19) 






Elimination rate constant of circulating neutrophils 


Calculated (Eq. 18) 




£ b 

max 


Maximum inhibitory effect of docetaxel on the proliferation of imma- 
ture cells 


83.9 




EC 50 (umol/l) b 


50% of maximum inhibition of docetaxel 


7.17 


77 


EC 50 (umol/l) (our data) 


50% of maximum inhibition of docetaxel 


28 c 


77 



Calculated: Values are calculated from other fixed parameters based on equations in Methods section. Simulated: Values are calculated by the in silico simulation 
based on the Monte Carlo simulation approach (see Methods section). 
a These PK parameters were taken from ref. 1 8. 
b These PD parameters were taken from ref. 17. 

c The EC 50 value was modified from the reported value (7.1 7) 17 to best explain our previous clinical results 6 because the very frequent MRP2 genotype was not 
taken into account. 

CV, coefficient of variation; PD, pharmacodynamic; PK, pharmacokinetic. 



/ = j60mg/h/m 2 (0~1h) ^iT^k = ^ x Transit, - /c (r x Transit 2 (15) 

[Omg/h/m 2 (1h~) dt 



^! = /c Prol x Prol x f^J x (1 - E Drug ) - k tr x Prol (10) dt 



dTransit 3 



: k tr x Transit 2 - /c tr x Transit 3 (1 6) 



, f dProl ^ 6t 

Prol tr [^T ) ( ] 



dCirc ^ x -fransit 3 - k c]rc x Circ (1 7) 



^=kcnc (18) 



Ed '»° = ec^Ci (12) MTT = ^T (19) 

X Detailed descriptions of all parameters are summarized 

C A = -j- (13) in Table 2. The means and coefficients of variation for each 

1 parameter are summarized in Table 2. 

To determine the prevalence of severe neutropenia in 

— ransi 1 = /c tr xProl-/c tr x Transit ( 14 ) specific populations by in silico simulation, each parameter 

d t in the PK/PD model for virtual subjects was generated as 
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described. 22 The distribution of all parameters was assumed 
to be lognormal. The number of virtual subjects generated 
in one simulation was set at 500. After entering the param- 
eters for each individual into the PK/PD model, the time 
course of neutrophil concentration in blood was simulated 
based on the differential equations described above using 
the Runge-Kutta-Gill method (Figure 4c). The calculation 
was performed with MULTI-RUNGE program coded by Visual 
Basic in Excel 2007 (Microsoft Corporation, Redmond, WA) 
as described. 22 

The grade of neutropenia was judged from the nadir based 
on the National Cancer Institute Common Toxicity Criteria 
version 2.0 (nonadverse drug reaction, >2.0x10 9 /l; grade 
1, 1. 5-2.0 x10 9 /l; grade 2, 1 .0-1 .5x 10 9 /l; grade 3, 0.5- 
1.0x1 0 9 /l; and grade 4, <0.5 x 1 0 9 /l). 

Our previous pharmacogenetic study showed that the 
intrinsic clearance of telmisartan in heterozygotes with the 
SLCOW3 rs1 1045585 mutation decreased to 70.5%. 19 
On the basis of this observation, we calculated the sys- 
temic clearance of docetaxel in homozygotes of SLCOW3 
rs1 1045585 using the following equation: 



CL tot - CL h 



Q h xf B x CL int h 



(20) 



where CL tot and CL h represent the systemic and hepatic clear- 
ances, respectively. Because docetaxel is eliminated pre- 
dominantly by hepatic metabolism, CL tot can be regarded as 
CL h . Q h represents the hepatic blood flow rate (20.7ml/min/ 
kg), f B represents the protein-unbound fraction of docetaxel 
in blood, and CL inth represents the intrinsic hepatic clearance. 
As for MRP2, because there is no information about the 
functional alteration in MRP2 caused by ABCC2 rs1 2762549 
mutation, we optimized the EC 50 values by minimizing the fol- 
lowing objective function (Obj). 



Obj = 



OR, 



wt/wt+wt/mut,obs 



-OR 



wt/wt+wt/mut,calc 



OR 



wt/wt+wt/mut,obs 



^wt/wt+mut/mut,obs " 



^wt/wt+mut/mut,calc 



OR 



wt/wt+mut/mut,obs 



OR 



wt/m ut+wt/m ut , o bs 



-OR 



wt/mut+wt/mut,calc 



OR 



wt/mut+wt/mut,obs 



(21) 



where OR obs and OR calc represent the clinically observed OR 
and calculated OR, respectively. 
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Study Highlights 



WHAT IS THE CURRENT KNOWLEDGE ON THE 
TOPIC? 

The risk of severe neutropenia caused by 
docetaxel was significantly associated with 
SNPs of SLCOW3 and ABCC2. Although 
docetaxel is a substrate of OATP1B3 and 
MRP2, the quantitative importance of their func- 
tional change in docetaxel-induced neutropenia 
is not well understood. 

WHAT QUESTION DID THIS STUDY ADDRESS? 

S The possible role of OATP1 B1 and MRP2 in the 
docetaxel-induced neutropenia was suggested, 
and its risk induced by SNPs of SLCOW3 and 
ABCC2 was quantified using PK/PD modeling. 

WHAT THIS STUDY ADDS TO OUR KNOWLEDGE 

y This study suggests that OATP1B3 determines 
the systemic exposure of docetaxel, while MRP2 
regulates its local concentration in the hemato- 
poietic cells. The ABCC2 rs1 2762549 is thought 
to decrease the MRP2 function to 32% to explain 
clinical observations, assuming that OATP1B1 
function decreased to 41% from our previous 
pharmacogenetic study of telmisartan. 

HOW THIS MIGHT CHANGE CLINICAL 
PHARMACOLOGY AND THERAPEUTICS 

y Docetaxel-induced neutropenia can be partly 
predicted from SNPs of SLCOW3 and ABCC2 
genes, and this will help in the careful manage- 
ment of patients with docetaxel-induced neutro- 
penia, depending on their genotype. 
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